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ABSTRACT: Individual-based models describe in- 
teractions between individuals in the population. Dy- 
namics of the number of individuals in the population 
is a byproduct of these interactions. Individual vari- 
ability plays an important role in this type of models. 
Individual-based models can be formulated in the form 
in which individual variability is replaced by average 
values. In contrast, classical models describe dynam- 
ics of population density directly and assume that all 
processes in the populations can be expressed as func- 
tions of the density. Numerical results of classical and 
individual-based models can be similar but interpreta- 
tions almost always are different. 


PE3IOME: VnansugzyanbHo-opHeHTupoBaHHble 
MOJEJIH OMMCbIBAaIOT B3AMMOJelCTBHA OCObe;i B NONY- 
JAHH. B sTHxX MOJeAX MHHAMHKA YHCJICHHOCTH OCO- 
Geğ B NONyYIAQHH ABJIAeTCA TOOOUHBIM IMpPOYKTOM 
ONHCÞIBACMÞIX B3AHMOMCĂCTBHĂ, a MHIMBHyaibHad 
H3MEHYHBOCTb HrpaeT BAXHYIHO pob. Mn quBuyyalip- 
HO-OPHeHTHpOBaHHBle MOEH MOXKHO nepeþopMyIH- 
POBaTb TaK, YTOŐbI MHAMBUyasIbHad U3MCHYMBOCTb 
Oba 3aMeHeHa CpemHHMH BEIHYHHAMH. Hanporus, 
KJIACCHYECKHE MOMEN HeEMOCpeACTBEHHO ONMCHIBAIOT 
JHHAMHKY YMCJICHHOCTH NONYyJIAQMĂ H MpHHUMAHT, 
YTO BCE MOMYIAWHOHHbIe NpPOHECCHI MOryT OBIT PyH- 
KUHel MWIOTHOCTH. UNCeHHbIe pe3yIbTaTbI OOOHX TH- 
IOB MOJeueH MOryT ObITb CXOHBIMH, HO MX HHTep- 
TipeTaljuu MOUTH BCerya pasyIM4HBl. 


1. Introduction 


Let us consider a population of a parthenogenetic 
insect species with a simple life cycle. Females lay 
eggs in autumn and die. Eggs overwinter in soil and 
hatch during the next spring. Larvae feed and grow. 
Adult females emerging in autumn do not feed, they 
only lay next generation of eggs and die. 


Even such a simple case can be complicated by 
some hidden assumptions. The first question is the 
survival of eggs, larvae and adult females. Not all the 
eggs laid in autumn survive until spring. Similarly, not 
all the larvae survive during summer and not all the 
females, which emerge from the pupae, survive until 
egg laying. However, to make the following consider- 
ations simpler, let us assume that all individuals, start- 
ing their lives as eggs, survive until egg laying time in 
autumn. This assumption enables us to leave the prob- 
lem of mortality aside and concentrate at this step on 
the reproduction process. We will come back, howev- 
er, to the mortality problem in the end of this paper. 

If we know the number of females in a given gener- 
ation and the number of eggs laid by each female, it is 
easy to calculate the number of females in the next 
generation. Let us assume that in generation ¢ there are 
N females. If i-th female lays r, eggs, the number of the 
females in the next [¢ +1] generation will be equal to 
the sum ofr, for all N females: 

Nf (1) 
where summation is fori=1,...,N. 

Two main questions arise now. The first one con- 
cerns the situation within the given generation. What 
are the values of r? Are they the same for each female 
or different? If they are different, can we predict their 
values? The second question extends for further gener- 
ations. If we want to calculate the number of females in 
subsequent generations, we have questions about their 
r values: are they the same or different as for the 
females in the present generation? If different, can we 
predict these values? 

This paper compares the possibilities of individual- 
based and classical modeling to describe dynamics of 
biological systems.It describes limitations and pitfalls 
of mathematical modeling in ecology. But it also stresses 
the usefulness of constant efforts to come with mathe- 
matical descriptions closer and closer to the reality. I 
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will consider an example of a simplest model of popu- 
lation growth — the exponential equation and its ex- 
tensions. 


2. Exponential equation — a short time scale 


Let us assume that all females are equal. They lay 
the same number of eggs r. Now, the equation (1) can 
be transformed to the following form: 

N,,=7N D 

However, a more interesting situation is when not 
all the females are equal. Let us first assume that in 
generation ¢ there are N' females, which lay r, eggs, 
and N? females, which lay r, eggs (of course N’ + N? = 
N ). The number of females in the generation [¢+1] will 
be described by the equation: 

N= tN Er N? (3) 

We can calculate r,,, — the mean (arithmetic) num- 
ber of eggs laid by a female from generation: 

Pa AN EN IN ENS) (4) 

The number of females in generation [ź +1] can be 
calculated using the mean number of eggs laid by a 
female from generation t: 

N =r aN, (5) 
Equation (5) i is in the form similar to equation (2). 

Nothing will change if we assume that all females 
are different and have different r, values. Once again 
we can calculate NV, using equation (1) or equation (5) 
and the result of calculations will be the same. Howev- 
er, in this later case we should use the following equa- 
tion for calculation the mean number of eggs laid by a 
m 

ma Èr) N, (6) 
d summation is for i = lysis N, 

There are two approaches to calculate the number 
of females N in the next [¢+1] generation. The first 
one (represented by equation (1)), which can be named 
individual-based, uses individual numbers of eggs laid 
by each female from generation ¢. The second ap- 
proach (represented by equation (5)), which can be 
named average-based, uses the average number of eggs 
laid by a female from generation ¢. Both approaches 
give the same result of N, and need the same amount 
of information to perform calculations. To use equa- 
tion (1), we need data on individual progeny produc- 
tion of each female. To use equation (5), we need 
average value of progeny production in a population. 
But to calculate the average number of eggs per fe- 
male, we also need individual information about the 
number of eggs laid by each female. Thus, both ap- 
proaches — individual-based and average-based — 
are entirely equivalent when applied at a short time 
scale (equal to one time step of a model, e.g.one year in 
our case population). 


3. Exponential equation — a longer time scale 


Let us now calculate the number of females for 
several consecutive years. We assume additionally, that 


all females, both in the given and future generations, 
are the same and lay the same number of eggs equal to 
r. In generation ¢ we have N females, the number of 
females n generations later will be: 

N „=N, (7) 

When in each ofn future generations the conditions 
are different, but it does not produce reproduction vari- 
ability among females within the generations, we have 
the a equation: 

Ns, (8) 
where multiplication i is for j = 1,...,n and r, is number 
of eggs laid by a female from generation J: 

Can we also apply in this case the approach which 
uses average value of the number of eggs per female? 
Let us calculate r „— a mean (harmonic) value of the 
number of eggs per female in consecutive generations: 


aye (9) 
where multiplication i is forj = 1,...,2. Then: 

Nein lm YN, (10) 
because 


Ca) =(=, (11) 

Equations (5) ‘and (10) describe the number of fe- 
males in consecutive generations using average values 
of eggs per female. Both equations are equivalent to 
the individual-based model (1). Equation (5) makes a 
prediction for one generation in advance, equation (10) 
for several subsequent generations. Prediction of equa- 
tion (5) is exact. Equation (10) predicts exact number 
of females in the generation [t+], however predictions 
for generations different than [¢+n] can not be made 
with the help of harmonic average for generation [t+n], 
but have to use harmonic averages calculated for prop- 
er generations. 

Results presented in this section can be easily gen- 
eralized to describe the case of a species with an indi- 
vidual variability in female reproduction in each gener- 
ation. In this case, to estimate the number of females in 
the future generations, we have to calculate arithmetic 
averages ofr values for each generation and harmonic 
average of these arithmetic averages. 

Again both approaches — individual-based and av- 
erage-based — are equivalent: they give the same re- 
sults and they need the same information for model 
construction. There is only one difference between them 
but its nature is rather numerical. Individual-based mod- 
els (1) use integers as input data and give results — 
numbers of females in the next generations — which 
are also integers. Average values calculated in aver- 
age-based models (5) or (10) are in general real num- 
bers and give numbers of females in the next genera- 
tions in the form of real numbers. However, average- 
based approach does not use this fact. This effect can 
be easily deleted by using function Trunc in computer 
simulations. 


4. Individual-based modeling 


Individual-based approach is represented by equa- 
tion (1). Let us replace r by r; ,to indicate that across 
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the different generations the number of eggs per female 
can be different. Then we have: 
(12) 


Nae 
where ‘summation is for all females present in genera- 
tion ti i= 1,...,,. 

The individual-based model does not consider the 
number of females as its main point. It concentrates 
rather on r; values and tries to explain, are they the 
same or individually vary among females. The number 
of females is a byproduct which we obtain only after 
the proper description of ecological processes hidden 
behind r/ values — growth, assimilation of resources, 
competition between individuals, unequal partitioning 
of resources as the result of competition, and costs of 
life too, which also can be different for different indi- 
viduals in the population. 

Individual-based model (12) (and its average-based 
counterparts (5) and (10)) can describe not only an 
exponential increase of the number of females in the 
population, but also population oscillations, its stabili- 
zation, decrease or even the extinction of the popula- 
tion. The dynamics of the system depends on a distri- 
bution of r values among females in each generation. 


5. Fitting to the empirical data 


Let us assume that we have the numbers of females 
in several subsequent generations and expect that the 
data follow an exponential increase of a population. 
However, we have no data on the number of eggs laid 
by females. Fitting the parameters of mathematical mod- 
el to the data is the most often action in such a case. We 
fit parameter r, of the equation 

Na =r N, (13) 
to the data we have. 

Parameter r,, which value has been obtained by 
fitting, is in general equal neither to r „ nor to r „ It 
gives the best fit to the empirical data. The numbers of 
females in subsequent generations calculated by means 
of equation (13) with fitted value of parameter r, give 
on average the best approximation of experimental data. 


6. Classical modeling 


Classical approach to model the dynamics of eco- 
logical systems assumes that it is possible to describe 
the number of females in subsequent generations by the 
equation: 

N,,=8N, (14) 
Equation (13) is similar to equations (5) and (10) 
using a procedure of averaging values of parameters 
characterizing individuals. Parameter g is in some sense 
related tor, or to 7 „ Its values can also be obtained 
by fitting to empirical data. Individual-based models 
use specific submodels to calculate r values which are 
different for different patterns of interactions between 
individuals in the populations. Classical models use a 
different approach. They try to be as general as it is 
possible, to include all possible forms of interactions 


between the individuals and describe them by some 
general functions. 

The main point of classical approach to population 
dynamics is an assumption concerning equation (14). It 
says that g is a function of N: 


g= g(N, ) (15) 
For instance, assuming that 
g(N,)=1+(1-N,/K) (16) 


where K is a constant parameter, we obtain the well 
known logistic equation, used to describe populations 
in which the number of individuals becomes constant 
after a sufficiently long time. 

Assumption (15) is of a statistical or empirical na- 
ture, because function (14) can be obtained by fitting to 
empirical data. It is also a very simplifying assumption, 
because it compresses complex ecological interactions 
(e.g., which influence female fecundity) in the form of 
a rather simple mathematical function. Only a very 
small set (two or three) of different functions is used to 
describe population dynamics of species with very dif- 
ferent ecology. 

This assumption has two important consequences. 
Applying equation (14) with assumption (15) shifts us 
from the world in which population number is ex- 
pressed by natural numbers into another world in which 
it is expressed by real numbers. In that case N is no 
more the number of females, but becomes a density of 
females in the population. This automatically limits the 
application of the model to populations which are well 
mixed in space and without spatial effects. Sedentary 
species can not be described by the model (14). 

Another consequence is very promising at a first 
glance. Classical models directly make use of the fact 
that population density is expressed by real numbers. 
Assuming that function (15) is a real function, we have 
an immediate access to all mathematic techniques which 
is used to analyze difference equations. In contrast, 
individual-based models are always simulation mod- 
els. This is a strong limitation, which can partly be 
overcome by proper simulation procedures (see the 
ODD protocols in Grimm et al., 2006, 2010, as an 
attempt to solve this problem). It is much more elegant 
to use methods available for analysis of difference 
equations when we are modeling with classical models. 
Analytical solutions are often possible in this case. 
However, the number of possible solutions is much 
more limited than in the case of individual-based mod- 
els. These are asymptotic solutions, cycles or deter- 
ministic chaos. Individual-based models produce a wider 
range of possible solutions. Working with classical 
models we devote much more time to stability analysis. 
From the point of view of classical modeling, a stabili- 
ty of ecological systems seems to be much more impor- 
tant than it is in individual-based models. This empha- 
sis on the stability of systems described by classical 
models is a direct result of mathematical methods used 
to model the dynamics of ecological systems. The im- 
portance of stability of ecological systems in nature is 
apparently much smaller (see Rhode 2005, for examples 
and arguments on so called nonequilibrium ecology). 
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Table 1. Data used to construct the individual-based model (number of females and r values). Harmonic average values of 
r (forth column) have been used to calculate results of the model (fifth column). 

Tadsuna 1. JJaHHble WIA MOCTpOCHHA MHIMBMyasIbHO-OpHeHTHpoBaHHOM MOJENN (4HCIIO CaMOK H 3HAYCHHA F). 
TapmMonuueckad cpeqHad 3HAYCHHĂ r (YeTBEPTad KOJIOHKA) HCMONb30BaHA JIA pacueTa pe3yJIbTATOB MOesMpoOBaHHA 


(MlocIeqHAA KOIOHKA). 






































Year Number of r values Harmonic Number of females, 
females, averages of calculated values 
empirical data r values 
1 100 2.0 2.0 100 
2 200 2.0 2.0 200 
3 400 2.0 2.0 400 
4 800 1.25 1.7783 800 
5 1000 0.8 1.5157 1000 
6 800 1.25 1.4422 800 
7 900 1.0 1.3687 900 
8 900 - - 900 














6. Conclusions 


Let us illustrate the arguments above by an exam- 
ple. Fig. | presents the numberof females in consecu- 
tive years. Numerical values of the number of females 
are also presented in the second column of Tab. 1. 
Other parts of Tab. 1 include all elements necessary to 
construct an individual-based model or its average- 
based counterpart. Values of r are the same for all 
females in each generation, but r values in different 
generations are different. Tab. 1 also presents harmon- 
ic averages calculated starting from year 1 for the 
subsequent years 2, 3, 4, etc., and the number of fe- 
males calculated according to equation (10). 

Tab. 2 illustrates the construction of a classical 
model based on the same empirical data (number of 
females and r values). The second and the third column 
are the same as in the previous table. According to the 
tules of classical modeling, we have to present r values 


Table 2. Data used to construct the classical model. The 
same number of females as in Tab. 1 was used. The func- 
tion of the form (17) has been fitted to describe a relation- 
ship between r values and the number of females. Equation 
(18) was then used to calculate densities of females in the 
population in the consecutive years. 

Tammya 2. JaHHble Wit MOCTpoeHHA KiaccH4eckoli 
MOJEJIH, YHCJIO CAMOK Takoe XKE, KAK B TaOm. 1. DyHKuMA 
(17) annpokcumyposBana MIA ONMCAHHA CBA3H MEKMY 
3HAYCHHAMH r H YHCIIOM CaMOK. Ypasnenne (18) 6suIo 
3aTeM HCHOJIb3O0BAHO JIA pacueTa IJIOTHOCTH CaMOK B 
TIOMYJIANH B PAY MOC OBATeIbHBIX JIET. 





























Year Number of r values Density of 
females, females, 
empirical data calculated values 
1 100 2.0 100.0 
2 200 2.0 216.2 
3 400 2.0 432.3 
4 800 1.25 733.6 
5 1000 0.8 935.5 
6 800 1.25 928.5 
7 900 1.0 930.6 
8 900 - 930.0 

















as related to the number of females and then fit a linear 
function to this relationship as it is most often done 
(see Fig. 2). The following function gives the best fit: 
g(N,) = 2.3022 — 0.0014 N, (17) 
After substitution of equation (17) to equation (14) 
we obtain the classical model of the form: 
N, = (2.3022 — 0.0014 N) N, (18) 
Forth column of Tab. 2 shows densities of females 
in the population calculated according to equation (18). 
Fig. 3 compares the results of individual-based and 
classical models constructed for the same data set. The 
difference between the results of both models is not big 
if we compare their ability to predict numerical results. 
Individual-based model gives the exact data. The clas- 
sical model only approximates the data, because it uses 
the function which is an approximation of a relation- 
ship between the population parameters. However, there 
exists a much deeper difference in the interpretation of 
these two models. Using function (18) we silently as- 
sume that independently of the population structure, at 


Fig. 1. Number of females in the consecutive years, as based 
on data used for construction of individual-based and classical 
models. 

Puc. 1. HACIO camoK B pAAy NOCIENOBATEJBHBIX JIET, Ha OCHO- 


Be JJ@HHBIX JIA MOCTPOCHHA MHAMBUyasIbHO-OpHeHTHpOBaHHOli H 
KylaccHueckoH MOJENEÑ. 
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Fig. 2. Number of eggs per female in relation to the number of 
females in the population. It is assumed that female fertility does 
not vary within, but varies across the generations. Linear relation- 
ship in the form of (17) was fitted to the empirical data. 

Puc. 2. Wucno su, OTIO2KeCHHBIX OJHOÑ CAMKOÑ, B 3ABHCHMOC- 
TH OT YACHA CaMOK B NMONYIANHH. [IpunnMaercs, YTO IJIOMOBH- 
TOCTb CAMOK HeH3MeHHa B MpeEgeJaX reHepallH, HO MeHAeTCA MEK LY 
reHepaluamu. JImHeiiHas 3aBHCHMOCTS B BHIE Þyukumu (17) 6bua 
alllIpOKCHMMpoBaHa K IMIMpW4eCKUM JaHHBIM. 





Number of females 











Fig. 3. Number of females in consecutive years: A comparison 
between results of an individual-based (solid circles) and classical 
(open circles) models. 

Puc. 3. ACIO caMoK B paly NoceqOBaTeIbHBIX JIET: cpaBHe- 
HHE MEK Ly pesyJIbTATAMH HH TMBULyaJIbHO-OpHeHTHpoBaHHoli (yep- 
HbIe CHMBOJIbI) H KIACCHYECKOÑ (Gebre CHMBOJIBI) MOJeJel. 


a given density all ecological processes in a population 
are fully and unambiguously explained by this density 
value. In contrast, during the construction of individu- 
al-based models we remember that for a given number 
of individuals in population there can be a variation in 
their demographic structure. For instance, the same 
number of individuals can have different weight distri- 
bution, different variability with respect to other im- 


portant individual features, or different proportion of 
females and males. Individual-based models always 
take such effects into account. 

Working with the classical model we can always 
look for singular points and perform stability analysis. 
Solving the equation 

1 = 2.3022 — 0.0014 N, (19) 
we find that N, = 930.1 is the singular point. It is also 
the stable point, because the absolute value of the 
derivative of function (18) with respect to N is less then 
1 for N = 930.1. 

Methods of individual-based modeling do not al- 
low to obtain results which are the product of mathe- 
matical technique used to construct and analyze the 
model. They do not make conclusions which have no 
roots in properties of the real system. After all, the data 
we have used in the example presented here do not 
indicate that there exists a stability of the population. 
Classical models directly use the fact that the popula- 
tion density is expressed by real numbers. This opens 
the opportunity to apply mathematical methods to per- 
form stability analysis of the model equation. But as a 
consequence, it leads to the results which are the prod- 
uct of mathematical methods and may be far from the 
reality. 

In the beginning of the paper we have assumed that 
there is no mortality in our case population. The intro- 
duction of mortality into the individual-based model 
does not change much the model and its interpreta- 
tions. The mortality of the progeny decreases r values. 
The mortality of adult females changes only the distri- 
bution of r values among the females. Similar numeri- 
cal changes will follow the introduction of mortality to 
the classical model. The analysis shows that relaxing 
the assumption about zero mortality does not change 
interpretations presented in this paper. 
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